Equivalence of channel-corrected T-matrix and anomalous propagator approach 



O 

(N 

G 
3 



K. Morawetz 1,2 

1 Miinster University of Applied Science, Stegerwaldstrasse 39, 48565 Steinfurt, Germany and 
2 International Institute of Physics (IIP), Universidade Federal do Rio grande do Norte - UFRN, Brazil 

Any many-body approximation corrected for unphysical repeated collisions in a given conden- 
sation channel is shown to provide the same set of equations as they appear by using anomalous 
propagators. The ad-hoc assumption in the latter theory about non-conservation of particle num- 
bers can be released. In this way the widespread used anomalous propagator approach is given 
another physical interpretation. A generalized Soven equation follows which improves a chosen ap- 
proximation in the same way as the coherent potential approximation (CPA) improves the averaged 
T-matrix for impurity scattering. 
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Superconducting and Bose-Einstein condensation phe- 
nomena belong to one of the most exciting macroscopic 
effects based on microscopic quantum physics. The theo- 
retical description of both phenomena is one of the major 
activity in theoretical many-body physics. Superconduc- 
tivity is based on the pairing of two fcrmions which form 
a condensate while bosons provides a one-particle Bose- 
Einstein condensate. Both phenomena are characterized 
by possessing a singular channel in which the effect ap- 
pears. 

The self-consistent multiple scattering T-matrix in a 
many-body surrounding diverges near the critical tem- 
perature of the onset of the symmetry-broken phase, may 
it be pairing condensation for fermions or Bose-Einstein 
condensation for bosons. Though describing correctly 
the onset of pairing, the T-matrix does not provide the 
gap equation. This changes if an asymmetric breaking of 
the self consistency in the T-matrix is used, such that one 
of the two internal propagators is used self-consistently 
and the other non-selfconsistently. Then the gap equa- 
tion appears as the pole of the T-matrix. This was first 
observed by Kadanoff and Martin and used later on 
and is called Prange paradox [f|. It has remained 
puzzling since a seemingly worse approximation leads 
to better results. Recently it turned out that the re- 
peated collisions Q with the same particle are responsi- 
ble for this artifact. Repeated collisions of two particles 
in the same state are unphysical since the particles move 
apart from each other after the collision. If these re- 
peated collisions with the same state are removed from 
the T-matrix the correct gap equation appears and the 
condensate can be described without asymmetrical ad- 
hoc assumptions about selfconsistency. The advantage 
of eliminating only the contributions of single channels 
as proposed in Refs. i and is that the formation of 
pairs and their condensation can be described within the 
same approximation. This has also resulted in the de- 
scription of different phases in interacting Bose systems 
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On the other hand there exist a well established theory 
to describe systems with condensates in terms of anoma- 



lous functions, for review see [8]. Let us consider bosonic 
particles which can form a condensate either bosons or 
paired fermions. If the system processes a macroscopic 
number No of such particles in the condensate repre- 
sented by channel i, the expectation value of the creation 
operator of that state is very huge and to a good ac- 
curacy the creation operators [ao, a q] « commute with 
each other as well as with all other states Q . The number 
of particles is considered as non-constant if the conden- 
sate is thought as a reservoir, since scattering off and on 
the condensate may create/destroy pairs. Therefore the 
anomalous Green function G\i = ^(Ta-k^k) are nonzero 
besides the normal Green function Gn — j(Taka~£) de- 
scribing the simultaneous excitation of a pair. This ad- 
hoc assumption leads then to the description of the con- 
densate and the gap equation for pairing. Please note 
that the nature of the condensate remains quite different 
whether it is composed of bosons or paired fermions Q . 

The question is now, how the two approaches above 
are related. In the first theory correcting the T-matrix 
we consider only microscopic correlations while the same 
result is obtained by the second approach where one as- 
sumes ad-hoc from the beginning anomalous functions. 
In the second approach the symmetry of the theory con- 
cerning particle conservation is obviously broken while in 
the first approach it remains conserving. So it seems to 
be worth to understand the relation between both ap- 
proaches. 

The aim of the present paper is to show that indeed the 
first approach leads to the same structure of equations 
like the second one yielding expressions for the anoma- 
lous functions without assuming them. Indeed it will be 
shown that the theory is somehow overdetermined by the 
second approach in that one can work with half the num- 
ber of equations if following the first procedure. In view 
of this, the virtue of the first approach consists in giving 
the anomalous function assumption a further microscopic 
meaning since they can be derived from theory. 

Let us therefore shortly sketch the structure of the first 
approach. We split the selfenergy into different chan- 
nels, E = Y) , Ej, where we assume for simplicity that we 
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have only one singular channel i where the condensate 
appears. The unphysical multiple scatterings with the 
same channel are concerning a single channel and vanish 
in the thermodynamical limit. Therefore this deficiency 
does not matter in normal matter. If we have a singular 
channel due to the condensate, however, this correction 
becomes essential. We have to subtract this process, i.e. 
we define the subtracted propagator 



G^ — G G^E^G 



(1) 



or G 1 = G^ 1 — Ej . Using the standard Dyson equation 
Gq 1 = G -1 + E we obtain the relation 



G\ — Go + Go(E — Ei)G^ 



(2) 



which shows that in this propagator the own selfenergy 
channel is subtracted, E' = E — Ej. Now we consider 
a general T-matrix which represents the selfenergy as 
£' = TjG where the channel T-matrix Tj as two- 

particle function is closed by an backward propagator G. 
In the singular channel we subtract the repeated inter- 
action within this channel. This is achieved by closing 
with the subtracted propagator Ej = TjG^. Now we can 
rewrite the Dyson equation as 



G-i fi— 1 
— o 



Go" 1 - £' 



= G 1 - E' - TiG\ 

= Gq 1 - E' — Tj (Go -1 - E' 



(3) 



where in the last step we have used ©. Finally we 
rewrite ([3]) to obtain the full propagator in momentum- 
energy (Matsubara) representation p — (p, w n ) 



G(p) = 



[G^{p) - £'(p)][Go x (p) - E'(p)] - T^p) 



(4) 



Remembering £' = E — TiGi leads immediately back to 
the Dyson equation G = Go/(l — EGo). Therefore it is 
an exact rewriting so far. 

Now we take into account the explicit form of the free 
propagator Gq 1 = uj — and Gq 1 = —lj — e_fc and call 
the "proper" selfenergy 



Eii(p)=E'(p). 



(5) 



Further, we observe that the T-matrix in the singular 
channel is separable flol - [l2| and can be written Tj(p) = 
=FA(p)A(p) for bosons/fermions respectively. Now we 
can define the "anomalous" selfenergy as 



Ei 2 (p) = A(p) 
such that the propagator (j4]) takes the form 
oj + e_ fe + En 



Gn 



(u + e-fe + Eu)(o; - e fe - En) ± Ej 



(6) 



(7) 



We add now an auxiliary quantity called "anomalous" 
Green function 



Gi 



T-Ei 



(u + e-k + Sii)(cj - e k - En) ± Ef 2 



(8) 



and can write in such a way the two equations in matrix 
form 



G = G° + G°EG 



(9) 



with 



G = 



Gn G12 
G12 ±Gn 



G = 



Go 
±Gi 



En E12 
E12 ±En 



(10) 



for bosons/fermions respectively. These are exactly the 
equations for anomalous propagators derived first by Be- 
liaev 13( for bosons. For fermions these are the Nambu- 
Gorkov equations 

In other words separating a singular channel from the 
selfenergy avoiding repeated collision within this channel 
leads immediately to propagators which have the Beliaev 
form for sosons or the Nambu-Gorkov form for fermions. 
We see that adding the auxiliary quantity ([8]) is not neces- 
sary. All information we have derived without this quan- 
tity and it was added here simply to show the same struc- 
ture of theory appears as provided by the approaches 
with anomalous functions. In this sense the theory of 
anomalous functions is overdetermined. We should note, 
however, that the anomalous propagator G12 describes 
the order parameter. This anomalous propagator ap- 
pears as a result of the theory here and is not assumed 
from the beginning as done usually. 

Now that we have clarified that the anomalous propa- 
gator is an exact rewriting of the Dyson equation if one 
correct a channel of selfenergy for repeated collisions we 
might ask what kind of equation such channel-corrected 
selfenergy obeys. This will lead us to a generalization of 
the Soven equation (IBj . The Soven equation was pro- 
posed to describe impurity scattering in terms of an ef- 
fective medium and resulted in the coherent potential ap- 
proximation (cpa) (lei-Hi 1 . This CPA improves the av- 
eraged T-matrix [19| with respect to better analytic prop- 
erties and a wider range of applications. It has turned out 
that the averaged T-matrix is the uncorrected channel 
while the CPA is equivalent to the channel-corrected ap- 
proximation. Here we will present the same idea of chan- 
nel correction but applied to the two-particle scattering. 
This will lead to a general Soven equation which allows 
to improve a chosen approximation scheme in the same 
manner as the CPA improves the averaged T-matrix ap- 
proximation for impurity scattering. 

Let s assume in general the defining equation for the 
channel T-matrix in terms of the potential V and a block 



T; 



V + VKjT, 



(11) 
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FIG. 1: Proof of relation (|16[) for singular channel T-matrices 
which become separable. 



covering both the singular channel j = i as well as the 
normal channels j ^ i. In general this equation is a two- 
particle one which is reduced to the one-particle selfen- 
ergy by closing the upper line with the backward propa- 
gating Green function £' = TjG and for the singular 

channel Ej = TiG^. In the following we will denote ex- 
plicitly by which function the upper line is closed. All 
other products are understood as operator products of 
one-particle functions. For the above mentioned aver- 
aged T-matrix approximation one has K = and a 
closing by G — > c in terms of the impurity concentration 
c. In two-particle ladder approximation one would have 
the form Ki = G ■ Gs. In the following we consider the 
general block K such that any more refined approxima- 
tion can be chosen. 

With ([T]) and Ej = TiG^ we can write 

1 = (1 - D)G- 1 G K + DG^G^ G^Gi^G^G (12) 

where we have added and subtracted an operator D 
which will be determined later by convenience. Now it 
easy to proof with the help of the separability of the 
singular channel Ti = =fAA that the following relation 
holds 



TiG$G — TiGG^. 
Indeed we have with the help of (UJ 

TiGG^G 1 — TiG — TiGG^TiG^ 

and 



TiG^ = TiG — TiG^iG. 
Since we have the identity 



TiGipiG — TiGG^TiG^ 



(13) 



(14) 



(15) 



(16) 



as shown in figure [TJ the equations (fl"4"|) and (|15l) are 
identical and relation (TT5]) is proved. 
Therefore we can write for (1121) 



Multiplying ([T]) from the right with TjG^ we can find 

G^G^G + E l G } 7K? = T~G = VG + VKTiG(18) 

where we have used (fTTj) for the second identity. For the 
last term we define now an effective potential Vi 



ViG^TiG = VKTiG (19) 
with the help of which we can invert (fT5|) 

G k T~G=(G- 1 + ^ l -V l )~ 1 VG. (20) 
Using this in (|17p we arrive at 
1 = (1- D)G- 1 G K 

^l + {^ l -V t )G]- 1 [(l + i: i G)D-V l GD + VGG}G- 1 G k 

(21) 

Now we choose the operator D such that the last two 
terms cancel each other, i.e. 



VGG = VjGD. 



(22) 



Using flU in the form G~ 1 G^ = (l+EjG)^ 1 and subtract- 
ing from HU the structure 1 = (1 - D)A~ 1 A + B~ l BD 
we obtain finally 

(l-DXl + EiG)"^ = 

[l + (E l -U)G]- 1 [^G^G- 1 -(l + E ! G)^(l + S l G)- 1 E J ]. 

(23) 



Together with the operator ([22]) and the effective po- 
tential Vi defined by (fl9|) this is the desired generalized 
Soven equation. It is written in operator form which be- 
comes an algebraic equation in the appropriate represen- 
tation. In the operator form it is even valid in nonequi- 
librium and its time ordering can be treated e.g. in the 
framework of generalized Kadanoff and Baym formalism 
[20| | with the help the Langreth/Wilkens rules 21 1. 

Though introduced merely for mathematical conve- 
nience, the operator D corresponds to the concentration 
for impurity scattering and the effective potential corre- 
sponds to the effective potential in CPA. Therefore D is 
called concentration operator hereafter. 

Let us illustrate this with the help of special cases. 
Choosing the averaged T-matrix approximation we have 
only one-particle functions, K — G^ and the closing by 
the concentration as a c-number G — ¥ c, such that we get 
from (HU) Vi = V which gives with (J22j> D = c such that 
the standard Soven equation [la ] appears 



(1-c) 



S<(p) 



V(p)-Zi(p) 



l + Ei(p)G(p) 



l + [Xi(p)-V(p)}G(p)- 



1 = (1 - D)G- 1 G i ^+ (DG- 1 + G~ 1 G(T i G)G{ (17) 



(24) 
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As a second, so far not known, example we give the 
explicit expressions for the two-particle T-matrix. Then 
K = G ■ is a product in spatial coordinates. Fourier- 
transform of the difference coordinates and gradient ex- 
pansion reveals then the structure [p = (p,cj p )] 



[1-D(p)} 



l + Ei(p)G(p) 



Dip) 



Vi(p)-Ei(p) 



l+[Ei(p)-Vi(p)]G(p) 
(25) 



The concentration operator takes the form 



D(p) 



K(P) 



Eg(p) 

Vi(p) 



(26) 



noting the Hartree selfenergy Ejj and the effective po- 
tential reads 



VKT t G TiG - VG 



1 



TjGGx 



Ti GGi 



g4p) 



Eg(p) 
Si(p)G(p)' 
(27) 



Here we have used ([TT|) for the first equality and (|T3l) for 
the second one. The channel T-matrix reads explicitly 



E,(p) =J> 4 



P-P P-P 



p+p,o; p + G\(-p) 



(28) 



It is instructive to see that introducing equations (|26p 
and ([27]) into (|25|) leads indeed to an identity. This is due 
to the fact that we have assumed that all quantities like 
T-matrix and selfenergy are known exactly. In case that 
we start with an approximation for the Green function 
or selfenergy we can use the above equation system and 
the generalized Soven form to iterate and to obtain ap- 
proximations for the channel corrected propagators and 
selfenergies. 

To summarize we have shown that the Dyson equation 
can be rewritten exactly into a form containing anoma- 
lous propagators and selfenergies if the repeated colli- 
sions in one channel are subtracted. The correction of 
the channel results into a selfenergy which obeys a gen- 
eralized Soven equation irrespective of the actually form 
of approximation used for the T-matrix. We have derived 
this identity here which allows to improve any chosen ap- 
proximation in the same way as the CPA does it with the 
averaged T-matrix approximation. The suggested proce- 
dure is to choose an approximation for the T-matrix and 



to calculate the self-consistent channel-corrected prop- 
agators by iteration of ([26]) and (12?]) and ([25]) . Due to 
the versatile appearance of pairing and condensation phe- 
nomena ranging from nuclear, to solid state up to plasma 
physics, the systematic improvement of any approxima- 
tion used and the underlying equivalences of the different 
approaches might be helpful. The actual application of 
the proposed scheme is reserved for further work. 
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